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Jets are studied in A-A collisions at RHIC and LHC with the goal to understand how they are 

H 

Cu affected by the medium and how they affect the medium. It is widely believed that hard- scattered 

partons lose energy when propagating through a medium before hadronizing. Partons losing 
enough energy may not even make it out of the medium as identifiable jets (although the mo- 
mentum will be shared among whatever particles are emitted). "Full" jet reconstruction attempts 
I— I to determine the partonic energy loss as well as possible changes in jet shape. Heavy ion col- 

lisions typically produce many unrelated particles within the jet "cone," and subtraction of this 
^—H background introduces significant uncertainties. A variety of techniques using high-pt particles, 

^ assumed to be leading particles from jet fragmentation, look for disappearance of jets and attenu- 

ation of jets relative to the reaction plane, as well as medium modifications such as Mach cones. 
Those techniques have considerable uncertainty due to subtraction of V2. In this paper we discuss 
minimum-bias jets observed at RHIC using two-particle correlations. We find that jets produced 
\0 in p-p collisions have interesting properties. Peripheral A-A collisions look like p-p collision. 
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ing until at a system-dependent centrality the number of particles associated with jets increases 



As we select more central collisions the number of jets increases following binary collision scal- 



substantially above this scaling. Near this transition centrality the jet aspect ratio — elongated 



transverse to the beam direction for low-energy jets produced in p-p collisions — becomes highly 
elongated along the beam direction in A-A collisions. 
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1. Introduction 

Our studies of two-particle correlations began in the late 1990's with the desire to understand 
heavy-ion collisions in an unbiased way [1]. At that time it was expected that central A- A col- 
lisions would be nearly thermalized with jets quenched and, if not completely absorbed, difficult 
to observe. Temperatures inferred from pt spectra could fluctuate event to event or from place to 
place within an event, the amplitude and volume of pt fluctuations depending on the nearness of 
the event to the quark-gluon plasma phase boundary. We wanted to observe these fluctuations in 
an unambiguous way. 

A number of techniques were considered, such as factorial moments [2] [3], wavelets [4] and 
entropy transport [5]. But those measurements are logarithmic in the scale, and we soon realized 
that the range of scales available to a detector like STAR was at best two orders of magnitude, 
making a measure linear in scale preferable. We eventually settled on measures related to variance 
differences and learned how to convert between correlations and fluctuations. [6] These correlation 
measures can be related to Pearson's correlation coefficient, making them more easily interpretable 
than fluctuations which are integrals over correlations and thus include correlation structure from 
many scales. 

There were some concerns about whether such a general approach could work, partly fueled 
by previous difficulty interpreting factorial moments [7] [8]. Also, it is easy for a detector artifact 
to produce a signal which may confuse the physics interpretation. With a variance-based measure it 
is possible to cancel detector artifacts in a natural way using mixed events. It was also believed that 
a generic two-particle correlation analysis wouldn't be sensitive enough to see correlations such as 
flow. Monte Carlo simulations of temperature fluctuations within nearly-thermalized events were 
performed and it was shown that our techniques could measure the fluctuations that were expected 
at that time. In fact we could distinguish the case where events were thermalized with each event 
having a different temperature from the case where events were nearly thermalized but having a 
temperature that varied with position [9]. When we analyzed collision data we saw considerable 
correlation structure, but no evidence for temperature fluctuations [10]. 

There are several structures apparent in two-particle correlations. These are typically different 
enough in angular size that we can disentangle them. There is a very sharp e^e~ peak due to 
y-induced pair production. The observed HBT peak is broader than the e^e~ peak but narrower 
than the same-side (SS) 2D peak due to minijets. When they are similar sizes we can use the 
fact that HBT is only observed for like-sign (LS) pairs while minijets have a strong unlike-sign 
(US) component. Back-to-back scattered partons are observed as an away-side (AS) ridge (as well 
as both partons contributing to the SS 2D peak) and this ridge is approximately independent of 
= ^1 — ^2 within the STAR TPC acceptance. For high pt particles we expect the AS ridge to be 
a Gaussian centered at tt on 0a = 0i — 02, but when that Gaussian is broad enough in 0a (as it is 
for most of our pt range) it has a cos (0a — ^) shape. We refer to cos (0a — tt) as an azimuth dipole. 
There is also a quadrupole observed, cos(20a). The dipole and quadrupole are orthogonal to each 
other on 0a and distinguished from HBT and the SS 2D Gaussian minijet peak by their long ranges 
in T]a. 

In the rest of this paper we emphasize the jet- like components of two-particle correlations. We 
start with a discussion of two-particle correlation measurement techniques, noting the equivalence 
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of fluctuations and correlations, and discuss how to interpret the multi-dimensional two-particle 
correlation space. Then we examine p-p collisions, the reference system for A-A collisions, and 
find interesting two-particle correlation structure. We next see how some of this structure is mod- 
ified in A-A collisions. Finally we make some comments about how pt correlations complement 
number correlations. 

2. Review of Methods 

Jets can be defined by a jet-finding algorithm which groups all particles that come from the 
fragmentation of a hard- scattered parton [1 1] [12] [13]. There are a few commonly used algorithms. 
It is an active area of research to determine which ones work best in a heavy ion collision envi- 
ronment. How well these algorithms work depends on the jet shape, which could depend on a 
jet-medium interaction as well as the fluctuating background of hadrons unrelated to the particular 
scattered parton. An alternative to explicit jet reconstruction is to study two-particle correlations. 
This approach makes no assumption about jet shapes and allows one to study lower-energy partons 
than explicit jet reconstruction algorithms permit. Indeed, below a few GeV, partons will only be 
able to fragment into two hadrons, and this is difficult to observe with jet reconstruction algorithms 
in heavy ion environments. Another advantage of two-particle correlations is that they allow study 
of inter-jet correlations as well as intra-jet correlations. 

The most common measure of the dependence between two statistical quantities is Pearson's 
correlation coefficient [14], the covariance divided by the geometric mean of variances, which is 
bounded by -1 and 1. We are interested in the structure of two-particle momentum- space correla- 
tions. We write the covariance Ap as a difference between an object and a reference, 

Ap {pi , ^2) = Psib {puP2)-pref {Pl , Pi) 

= P(2) {Pl , P2) - P(l) (^1 )P(1) {P2) 

where pp) is a two-particle density and P(i) is a one-particle density. In practice we bin quanti- 
ties, storing them in histograms. For the case of number correlations the bin (a,/?), representing 
correlations between particles at positions a and b, can be written as 

£a£bAp{n) (riant) {fia) (nt) {n-n)a{n-n)t 

where £ are the bin widths. For Pearson's correlation coefficient we divide by the geometric mean 
of the variances, Ap / -^cj^^cj^^ ^ Ap/ y/W^ = ^P/ ^JPref- We have replaced the variances in the 
denominator with the Poisson expectations so that all the physics is in the numerator. In a detector 
one must deal with efficiencies and acceptances. The quantity we actually use is thus 

We normally refer to this simply as Ap / ^ypref' We define p'ref = ^^^^^ ^P / P^^f 

the virtue of canceling efficiencies and acceptances. Besides being closely related to a standard 
correlation measure the quantity Ap/y/p^ is the correlation strength per final-state particle. If 
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Figure 1: Sequence going from measurement of fluctuations at STAR fuU scale to scaled fluctuations prob- 
ing dependence on size and finally to correlations inferred from fluctuations. There is an interesting structure 
in the correlations. 

A-A collisions are simply superpositions of N-N collisions it will be independent of centrality and 
collision system. 

Correlations and fluctuations are closely related, fluctuations resulting from an integral over 
the correlation structure. For the angular space the integral equation relating the fluctuations and 
correlations is [6] 

r5ri r5(j) 

Acj;^.^(5r],50)=4 / driA / J0a^(577,50;t7a,0a)^=(t?a,0a), 

^ ^0 Jo ^/pref 

where A<7^^.^(5t] , 50) is a variance difference measured in a "detector" of size (Sr/ , 50), the kernel 
K is an exactly known geometric factor and Ap / ^ypref is a two-particle autocorrelation. This form 
of integral equation is known as a Volterra equation. Writing this equation in terms of binned 
quantities we can evaluate the kernel explicitly (Eq. 5 of [16]) 

Here we use A instead of p to indicate it is a binned quantity. We see that in general fluctuations 
depend on the domain scale over which they are measured. When we measure fluctuations at a 
particular scale we have an integral of the correlations up to that scale. The difference in fluctu- 
ations between two scales is a measure of the integral of the correlations between those scales. 
We can measure the scale dependence of Acr^.^ on (5t],50) and invert the integral equation to 
infer -^=(t]a, 0a)- Inverting Volterra equations (to solve for correlations from fluctuations) is an 
example of an ill-posed problem. One must impose a regularization scheme to obtain a useful 
solution, analogous to applying a low-pass filter. There is a well-developed mathematical frame- 
work for this [15] but one always has to assess how much the regularization scheme affects the 
extracted correlations. An example of inverting scaled fluctuations to get correlations is shown in 
Fig. 1 [17] [18]. We note that one motivation for measuring moan-pt fluctuations was to determine 
event-to-event temperature fluctuations. The correlation structure we observe is inconsistent with 
temperature fluctuations, instead it reveals jet correlations. 

Fluctuations are less numerically expensive to compute than correlations but they are difficult 
to interpret. We can invert the integral equation to get the more-easily interpretable correlation 
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(a) (^A 5 0a) projection (b) (j^i,_y^2) projection (c) US, AS projection (d) US, SS projection 



Figure 2: Projections of the six dimensional correlation space onto two dimensional subspaces for p-p 
collisions. Panel 2(a) is the projection onto the angular subspace, (t7a,0a), with a "hard" pt cut imposed. 
We see an intra-jet correlation on the same-side (SS) and an inter-jet ridge on the away-side (AS); the ridge is 
due to the fact that the parton center of mass is not at rest in the laboratory frame. Panel 2(b) is the full space 
projected onto the (yti^yti) subspace. Panels 2(c) and 2(d) are projections of the AS and SS correlations of 
unlike-sign (US) pairs onto the (yti^yti) subspace. 

structure, but the issue with regularization will always be a source of uncertainty. There are also 
technical issues such as two-track resolution that are difficult to correct in a fluctuation measure. 
Having established a direct connection [6] we prefer to measure angular correlations directly. 

Each particle has three momentum components, so the two-particle correlation space is six 
dimensional. In a detector with complete azimuthal coverage the absolute position of the pair (0^ = 
01 + 02) does not matter. Only the difference 0a niatters. Within the T] acceptance of the STAR 
TPC it is approximately true that pair correlations are independent of T]^ = t]i + T]2, depending only 
on difference T]a. Averaging over 0£ and 7]^ reduces the dimensionality of the correlation space to 
four with no loss of information. 

We are still left with a four-dimensional space, two angle differences and two momentum 
magnitudes. It is common to use pt for the momentum magnitude, but since the yield falls steeply 
with Pt most of the pairs occupy a small corner of the relevant {pn^Pti) space. We use yt = 
In [{mt + Pt) /mj^] for better visual access to the low-;?^ correlation structure. Four dimensions is 
impossible to visualize so we project onto the 2D angular subspace, (t/a, 0a), which has the form 
of a joint autocorrelation, or the 2D {yti^jti) subspace. It is instructive to make cuts on the space 
that is projected out. In Fig. 2 we show examples of projections onto the (t]a,0a) and (^^1,3^^2) 
subspaces. For the (^^1,3^^2) subspace we show the projection over the full angular space as well as 
for the AS and the SS. 

3. Two-particle correlations from proton-proton collisions 

In this section we examine the two-particle correlation structure of particles produced in p-p 
collisions. First we make a connection with a spectrum analysis in which events with different 
numbers of observed charged particles {rich) were analyzed separately. It was observed that one 
can define a "soft" component. When the soft component is subtracted from each of the rich spectra 
what remains has a shape and position independent of rich, but an amplitude proportional to rich — 
the so-called "hard" component [19]. When plotted on yt the hard-component shape is Gaussian 
with a peak at yt = 2.7 or pt ^ I GeV/c. This is the same position of the peak we find in the 
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(a) SS, LS (b) SS, US (c) AS, LS (d) AS, US 



Figure 3: Projections of two-particle correlations onto (yti^yti) with cuts on angular space. The AS parton 
fragmentation region is independent of charge sign, the US pairs have additional soft component. The SS 
parton fragmentation region is restricted to US pairs. 

projection of the two-particle correlations onto (^^,^^2) as shown in Fig. 2(b). In Fig. 5 we will 
see that for low-ytz = yti -\-yt2 pairs the angular correlations are dominated by soft physics while 
for high-ytz pairs the angular correlations are dominated by hard scattering. 

We examine (^^1,3^^2) in more detail by cutting on SS and AS pairs as well as looking at LS 
and US pairs. These four combinations are shown in Fig. 3. We refer to the region around yt ^3 
(pt ^ 1.5 GeV/c) as the fragmentation region. We see that for the AS this region is independent 
of the charge combination. For SS pairs the fragmentation region is restricted to US pairs. These 
correlations are dominated by low-Q^ ^ 3 GeV partons which primarily fragment into hadron pairs, 
the pair being charge neutral. 

We have already seen in Fig. 2(a) that when projecting from the fragmentation region of 
{yn^yti) the angular correlations have clear hard scattering (jet) structure. We can also ask how 
the shape of the {yti.yti) correlations should look if it is really due to fragmentation. One can 
parameterize parton fragmentation functions (FFs) in a universal form [21]. This allows us to con- 
struct a 2D plot with the parton momentum along one axis and the fragment momenta along the 
other axis as shown in Fig. 4(a). A slice parallel to the fragment momentum axis is a fragmentation 
function. The amplitude of these slices is determined by the underlying parton spectrum. For large 
parton momentum the mode of the fragmentation function is well determined by pQCD [20] [21]. 
As the parton energy decreases the number of hadrons it fragments to decreases until for some 
energy 3GeV) there are only two fragments (a parton may turn into a single hadron "fragment", 
but we are studying two-particle correlations). We do not observe the parton, only the hadrons. So 
we symmetrize the parton versus fragment plot to get a fragment versus fragment plot, integrat- 
ing over parton momentum as shown in Fig. 4(b). At low ytz the yield is dominated by partons 
fragmenting to two hadrons and is peaked along the diagonal. As the fragment yt increases the 
importance of parton fragmentation to three, four and more particles increases. The maximum in 
the (yti^yti) plot deviates from the diagonal and follows the line of modes predicted by pQCD. 
This expectation is sketched in Fig. 4(b). Compare this to the measured (^^1,3^^2) correlation for SS 
fragments in Fig. 4(c). 

We can also go in the other direction, from the data to the line of modes. We take the (yti^yti) 
plot and project the conditional slices shown in Fig. 4(c) onto yt. For relatively large yt these are 
well described by a Gaussian with the same width determined from the spectrum analysis [19]. 
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(a) parton-fragment (b) fragment-fragment (c) conditional cuts (d) conditional slices and locus of modes 



Figure 4: Panel 4(a) shows the parton momentum versus fragment momentum. A vertical slice projected 
onto fragment momentum is a fragmentation function. In panel 4(b) we symmetrize, integrating over parton 
momenta to get a figure we can compare to measurement. We note that at low momentum the fragment yield 
is dominated by partons fragmenting to two hadrons and is peaked along the diagonal. For larger momentum 
we can have parton to three or more fragments and we see the mode following the pQCD expectation. 
Panel 4(c) is the measured SS (yn^yti) correlation which should be compared with 4(b). In panel 4(d) we 
show conditional slices for the cuts defined in panel 4(c), fitting them to get the modes of the fragmentation 
functions. These closely follow the expectation sketched in panel 4(b). 

Using this width to fit all yt slices works well down to yt of about 2.5 (pt ^ 1 GeV/c). A plot of the 
Gaussian centroid versus yt slice starts off along the diagonal then curves, following the expected 
pQCD line of modes. This is shown in Fig. 4(d) which should be compared with Fig. 4(b). The 
interpretation that the area around yt ^3 is due to fragmentation is supported not only by the shape 
of the angular correlations but also by details of measured fragmentation functions. 

We now turn to the angular correlation dependence on the (^^1,3^^2) space. In Fig. 5 we show 
the angular correlations for LS and US pairs for "soft" (ytz < 2.7) and "hard" (yti^yti > 2) pairs. 
We see a pronounced HBT signal in the LS pairs at (0,0), especially for the soft cut. The largest 
component of the soft US pairs is a ID Gaussian on T]a due to charge ordering in projectile nucleon 
fragmentation along the beam axis. Imposed on this is a dip centered at (0,0) due to momentum 
conservation (when there are only a few particles in an event they will be biased against having 
the same direction, unless they are fragments from a single hard- scattered parton) and a narrow 
peak at (0, 0) due to e^e~ from 7-induced pair production. The hard pairs have an AS ridge that is 
independent of charge, the SS peak is primarily in US pairs. 

In Fig. 6 we examine this jet structure in more detail by making finer cuts on 3;^^. As ytz is 
increased the SS becomes sharper and the AS ridge becomes narrower. The width of the AS ridge 
is affected by a measure of the initial-state parton transverse momenta. The widths of the SS 
2D Gaussian are related to jV?] and jt^, measures of the transverse momentum with respect to the 
thrust axis of the fragmentation process [22] [23]. The relations between the widths and the jtS 
and kt is usually derived for the asynmietric case of a leading particle [24] [25]. In this analysis 
we use a formulation symmetric in the two particles [26]. In addition, since the most-probable 
parton has Q^3 GeV and the observed kt is typically around 1 GeV/c (with the jtS being not too 
much smaller) we cannot use small-angle approximations. When ytz is comparable to jt or kt the 
fragment distributions will be affected by kinematic limits. The results of fits to the SS peak are 
shown in Fig. 6. Panel 6(a) shows the widths in the rj/s, and 0a directions. Since yt is the log of pt 
the X-axis is essentially l/pt, the pt of the pair decreasing from left to right. We see that for this 
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(a) "soft" LS (b) "soft" US (c) "hard" LS (d) "hard" US 



Figure 5: Projections of two-particle correlations onto (r7A, 0a) space with cuts on {yti^yti) space. Soft pairs 
have yti -\-yt2 < 2.7 and hard pairs yti^yti > 2. For the soft pairs the LS has a strong HBT component while 
the US is dominated by local measure conservation such as momentum and charge. For the hard pairs the 
AS ridge is charge independent while the SS jet peak is primarily US. 

variable the T] and widths of the SS peak are approximately linear, but different from each other 
and dependent on charge combination. Panel 6(c) shows the inferred j^s, this time plotted on y^z- 
The difference in the jtri and jt(i) is greatest at low ytz- As ytz is increased they both increase until 
at large enough momentum they approach the observed perturbative j^-scaling value [23]. 

The width of the AS ridge can be used to infer kt, which has components from intrinsic initial- 
state parton momenta and initial- state radiation [27]. The momentum transfer between the scattered 
partons defines a direction q. When the initial kt of both partons is perpendicular to q but parallel 
to each other the AS ridge width is widest. The effects of kt also show up in the AS {yti^yti) 
correlations. When the ^^s of the two partons are parallel to each other and parallel to q the yt 
values for the fragments of one of the partons will be boosted while the yt for fragments of the 
other parton are reduced. This populates the off-diagonal region of the (^^1,3^^2) correlation. We 
show the AS (yti^yti) correlation in panel 7(a). In contrast, the SS correlation, shown in panel 7(b), 
is unaffected by kt and is considerably narrower. Typically, measured kt is about 1 GeV/c, nearly 
the same momentum as the minijet partons we are probing. It may be possible in some cases for 
the kt to boost the scattered partons enough so that they both actually emerge on the same side. 

We have seen how we can use standard correlation measures to extract information in p-p col- 
lisions. We find clear signatures of hard parton scattering at surprisingly low parton momenta. We 
have also seen that there is a nice correspondence between two-particle correlations and a spec- 
trum analysis which indicates a hard component appearing in perhaps 1% of NSD p-p collisions 
at 200 GeV [19]. We want to use this information as the baseline when analyzing A- A collisions, 
keeping in mind that the fraction of particles produced by hard scattering depends on centrality. So 
the superposition of N-N collisions should be weighted appropriately to account for the hard p-p 
spectrum component. 

4. Two-particle correlations in A-A collisions 

In this paper we only discuss angular correlations from symmetric A-A collisions. We have 
analyzed Au-Au collisions at 200 GeV and 62 GeV as well as Cu-Cu collisions at the same ener- 
gies. We see evolution with centrality from a correlation structure consistent with p-p collisions for 
the more peripheral events to a structure dominated by the SS 2D peak being strongly elongated 
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(c) jtj^ and jt(j) (d) SS peak for ytz ^ 6.5 



Figure 6: Panel 6(a) shows the fitted widths of the SS peak, the x-axis being essentially l/pt of the pair. 
The widths depend on the direction (rj or 0) and the charge combination but are approximately linear in this 
variable. Panel 6(c) shows the extracted jt values for the rf and (j) directions, having very different values for 
low yt but approaching the jt scaling value at high yt. Panels 6(b) and 6(d) are projections onto (t7a, (/>a) for 
low and high ytz cuts and plotted in 1 : 1 aspect ratio and graphically showing the change in the shape of the 
SS peak as a function of pair momentum. 

along T]a and an AS ridge consistent with back-to-back hard parton scattering. For mid-central 
collisions there is in addition a cos (20a) quadrupole component. Examples for a few centralities 
of 200 GeV Au-Au collisions are shown in Fig. 8. The other collision systems are qualitatively 
similar. 

We quantitatively describe the 2D correlation structure with a SS 2D Gaussian (with different 
T]a and 0a widths), a cos(0a — ^) dipole describing the AS ridge and a ID Gaussian on 7]^ [28]. To 
complete the description of p-p collisions we include a narrow exponential at the origin to describe 
e^e~ pairs from y conversions as well as an overall offset. We don't infer any physics from the 
exponential peak or the offset. This model function with ten parameters works well for the most 
peripheral and fairly well for the most central data. For intermediate centralities a cos(20a) term is 
required as well. 

Fit parameters are shown in Fig. 9. The most striking feature of the fit parameters is that the SS 
Gaussian amplitude closely follows the expectation for binary scaling up to some centrality then 
greatly exceeds it. The AS ridge described by cos(0a — Tt) tracks the amplitude of the SS Gaussian. 
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(a) AS (ytuyti) (b) SS (ytuyti) (c) "hard" (t7a, 0a) 



Figure 7: Effect of kt on (yti^yti) correlations. Panel 7(a) shows the AS (yti^yti), broadened when the 
kt of the two scattered partons are parallel to each other and the direction of the momentum transfer. The 
SS correlation, shown in panel 7(b), is not affected by kt. Panel 7(c) indicates how one can extract kt by 
measuring the width of the AS ridge. 




(a) 84-93% (b) 55-64% (c) 46-55% (d) 5-9% 



Figure 8: Samples of (77a, 0a) correlations for 200 GeV Au-Au collisions. Panel 8(a) is most peripheral 
(84-94%) while panel 8(b) is nearly central. The middle two panels are on either side of the sharp transition, 
55-64% and 46-55% centralities. 

The r/A width of the SS peak also increases greatly close to the same centrality. In contrast to the t/a 
width, the 0a width starts at the same value as we find for p-p collisions but then actually decreases. 
In contrast to the sharp transition of the jet amplitude the cos(20a) component has an interesting 
but unrelated evolution on centrality. 

We emphasize the deformation of the SS Gaussian peak by plotting the data after subtract- 
ing the multipoles determined by the fit (we also subtract the sharp exponential peak) in Fig. 10, 
plotting the data in a one-to-one aspect ratio. The most peripheral bin is actually elongated along 
the 0A direction. But as we increase centrality the peak becomes symmetric, then within a narrow 
centrality range becomes greatly elongated along t]a. 

In Fig. 11 we present an isometric view of the same data. We subtract the fitted dipole, 
quadrupole and sharp exponential peak from the data. The sharp peak affects only a few bins 
around the origin. Our parameterization of the data does a good job describing the 0a structure 
of the AS peak. There is a small but statistically significant T]a structure remaining on the AS. 
We note that the 0a dependence of the SS 2D peak plotted here can be decomposed into Fourier 
components on azimuth and would contribute significantly to an inferred quadrupole component. 
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Figure 9: Values of fit parameters for Au-Au colHsion systems. For the more peripheral data the SS peak 
amplitude follows the binary scaling expectation, exceeding it dramatically at a transition centrality that 
depends on the energy/collision system. The AS ridge follows this trend very closely and the 77a width 
of the SS peak also increases substantially at a similar centrality. In contrast the 0a width of the SS peak 
decreases with increasing centrality and the cos (20a) has an interesting but smooth evolution, already 
having a significant value before the transition centrality. 



Multipole subtracted, 84-93% centrality | 



Multipole subtracted, 55-64% centrality 



Multipole subtracted, 46-55% centrality | 



Multipole subtracted, 5-9% centrality | 







(a) 84-93% 



(b) 55-64% 



(c) 46-55% 



(d) 5-9% 



Figure 10: Samples of (r]A, 0a) correlations for 200 GeV Au-Au collisions after multipoles and sharp peak 
have been subtracted and plotted in a 1:1 aspect ratio. Panel 10(a) is most peripheral (84-94%) while 
panel 10(b) is nearly central. The middle two panels are on either side of the sharp transition, 55-64% 
and 46-55% centralities. 
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I Multipole subtracted, 84-93% centralily | Multipole subtracted, 55-64% centrality | Multipole subtracted, 46-55% centrality | Multipole subtracted, 5-9% oentrality 




(a) 84-93% (b) 55-64% (c) 46-55% (d) 5-9% 

Figure 11: The same centralities as Fig. 10 but plotted in a view similar to Fig. 8. We have subtracted the 
fitted dipole, quadrupole and sharp exponential peak from the data and rotated slightly. We see the dipole 
and quadrupole terms have exhausted the 0a dependent structure on the AS. 




(a) 84-93% (b) 55-64% (c) 46-55% (d) 5-9% 

Figure 12: Samples of (t] a, 0a) Pt correlations for 200 GeV Au-Au collisions. These are the same centrali- 
ties shown in Fig. 8. 



We believe this peak is due to hard scattering (it certainly is for peripheral collisions) and in any 
case this shape cannot be due to a medium response from a pressure gradient. 

5. Pt correlations 

We now return to the topic of pt correlations which we mentioned briefly in section 1 . For 
number correlations the Pearson's correlation coefficient quantifies how the number density at two 
locations vary with respect to each other. One could also ask how the total momentum densities 
are correlated. This turns out to be still dominated by number correlations. We instead ask how the 
mean ptS are correlated. Specifically, 

Ap {Pt-npt)a{pt-npt)h 

y/Pref V^a^b 

where pt is the mean pf of the parent population. We show examples of these correlations in 
Fig. 12. There is a qualitative similarity with the number correlations but there are quantitative 
differences. Even the most peripheral bin shows no indication of a ID Gaussian on T]a. The sharp 
peak at (0,0) is dominated by HBT, opposed to number correlations where e^e~ is a significant 
contribution. The broader peak centered at (0,0) is a 2D Gaussian for number correlations but has 
a catenary shape on T]a for pt correlations. 
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(a) 80-90% (b) 45-55% (c) 20-30% (d) 0-5% 



Figure 13: pt correlations for 200 GeV AuAu collisions for various centralities. These correlations have 
been inferred from the inversion of scaled fluctuations and we have subtracted the multipoles. The shapes 
of the peaks at (0,0) are different than those of number correlations and there is a "recoil" hole around the 
peak. 

We previously looked at mom-pt correlations when we were measuring pt fluctuations and 
found we could invert fluctuation scale dependence to obtain detector-acceptance-independent cor- 
relations [18]. We found that the SS 2D peak amplitude was nearly proportional to the mean 
participant path length (v) from the most peripheral to about 30% most-central, falling for the most 
central. As the centrality selection increased the SS 2D T]a width increased by about 60% while 
the 0A width decreased by about 30%. These trends are qualitatively similar to number correlations 
but quantitatively very different. In addition, there appears to be a "recoil hole" around the SS 
peak [18] which is not observed in number correlations. We also found that although Hijing does 
produce a SS 2D peak and an AS ridge in pt correlations these are nearly centrality independent as 
well as having the wrong detailed shape, with no indication of a recoil hole [18]. We are presently 
working on a detailed description of these correlations to characterize our direct pt correlation 
measurements. 

6. Summary 

We have presented a detailed differential analysis of minimum-bias jet systematics in p-p 
collisions. We have described low-Q^ jets using fragmentation function systematics. The widths 
of the same-side (SS) 2D peak depend on the ytz of the pair, changing from elongation along 0a 
to symmetry in (t]a,0a), eventually reaching a width described by kt scaling as ytz is increased. 
We found that kt not only affects the away-side (AS) angular correlations but also broadens the AS 
(ythyti) correlations. 

For A-A collisions we saw that the SS 2D peak elongates along T]a, the width increasing 
rapidly at a particular centrality (sharp transition). The amplitude of the SS peak is consistent with 
binary scaling for peripheral collisions and greatly exceeds binary-collision scaling starting near 
the same centrality where the T]a width elongation starts. The AS dipole amplitude closely follows 
the SS 2D Gaussian peak amplitude. In contrast, the 0a width of the SS Gaussian decreases slightly 
with increasing centrality. The quadrupole component is small for peripheral and central collisions 
but has already become significant before the transition centrality, where particle densities are low. 
There is no evidence for an opaque core. All partons, even low-2^ partons, are accounted for in the 
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final state. We are currently studying pt correlations which provide complementary information on 
jets. 
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